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Abstract 

The fuzzy disc is a matrix approximation of the functions on a disc which preserves 
rotational symmetry. In this paper we introduce a basis for the algebra of functions on 
the fuzzy disc in terms of the eigenfunctions of a properly defined fuzzy Laplacian. In the 
commutative limit they tend to the eigenfunctions of the ordinary Laplacian on the disc, i.e. 
Bessel functions of the first kind, thus deserving the name of fuzzy Bessel functions. 



1 Introduction 



Noncommutative geometry [1, 2, 3, 4] has suggested the introduction of fuzzy spaces, namely the 
approximation of the abelian algebra of functions on an ordinary space, with a finite rank matrix 
algebra, which preserves the symmetries of the original space, at the price of noncommutativ- 
ity. The first formalization of the idea of a fuzzy space was introduced by Madore [5], for the 
sphere. In his approach, a fuzzy sphere is a sequence of nonabelian algebras, generated (loosely 
speaking) by three "noncommutative coordinates" which satisfy XiXi = 1 and [xj,Xj] = KEijkXk, 
with k some noncommutativity parameter. The connection with the finite dimensional repre- 
sentations of SU (2) is immediate, and with k proportional to the inverse of the rank of the 
representation, one can think of recovering functions on the sphere in some limit. This limit has 
been made rigorous by Rieffel [6] at the level of topological "quantum spaces" , but here we are 
more interested in the approximations that the fuzzy matrix geometry can give to the space of 
functions (and hence of fields) defined on some space. This is efficiently done defining "fuzzy 
spherical harmonics" [7, 8], a basis for the matrix algebra, which in the limit converges to the 
ordinary spherical harmonics. Field theory on the fuzzy sphere is a very active field of research. 
A partial list of references is [9]. 

After the fuzzy sphere, similar approximations for higher dimensional spheres [10], for the 
torus, based on the noncommutative torus algebra [11, 12], and for projective spaces [13], have 
also been studied. In previous work [14, 15], we proposed a fuzzy approximations for the algebra 
of functions on a disc (see also [16, 17]). It is the first example of a fuzzy approximation for a 
space that classically has a boundary. The space is defined through a suitable truncation of the 
noncommutative plane. In this paper we examine the spectral properties of the fuzzy Laplacian, 
finding its eigenvalues and eigenvectors, and comparing them with the exact case and give a 
direct way to associate, to functions on the disc, their fuzzy counterpart, i.e. to approximate 
the abelian algebra of functions by a sequence of finite rank matrices. We introduce a basis 
for the algebra of functions on the fuzzy disc in terms of symbols of the eigenmatrices of the 
fuzzy approximation to the Laplacian. These symbols are seen to converge, as the rank of the 
matrices increases, towards a set of functions on the plane which coincide with the ordinary 
Bessel functions of first kind (eigenfunctions of the ordinary Laplacian with Dirichlet boundary 
conditions) for points on a plane inside a disc, and to zero for points outside this disc. These 
are given the name of fuzzy Bessel functions. 

The paper is organised as follows. In the next section we review the fuzzy sphere from a 
perspective which will render easy the generalization to the disc case. The main tool we use 
is the Weyl-Wigner formalism, and the concept of generalized coherent sates, which we briefly 
review in appendices A and B respectively. In section 3 we review the fuzzy disc. A Weyl- 
Wigner isomorphism is introduced in terms of functions on the plane, where noncommutativity 
is represented by a parameter 9. In this formalization, there is no natural concept of a sequence 
of finite dimensional Hilbert spaces, or finite rank matrix algebras, therefore a truncation in 
the algebra of operators is introduced, with respect to a specific basis in the Hilbert space. 
Constraining the dimension N of the truncation to the noncommutativity parameter, N6 = R 2 , 
one obtains then a sequence of finite rank matrices, converging towards an abelian algebra of 
operators, that approximates functions whose support is concentrated on a disc of radius R. 
On this sequence of states, fuzzy derivatives and a fuzzy Laplacian are defined and the spectral 
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properties of the latter are studied. Then in section 4 we study eigenvalues and eigenvectors 
and show that they (rather quickly) converge to the ordinary Bessel functions. 



2 The fuzzy sphere in the Weyl-Wigner formalism 

In this section we present the fuzzy sphere in the general context of the Weyl-Wigner formalism 
[19, 20, 21]. We first set up an isomorphism between a space of operators and a space of functions 
on a sphere. Since a sphere is the coadjoint orbit of the group SU (2), the basic tool to introduce 
this map is a system of coherent states, specialising the general arguments of appendix B. 

To define a system of coherent states consider the unitary irreducible representations for 
the group SU (2). On each finite dimensional Hilbert space C^, with N = 2L + 1 -here 
(L = 0, 1/2, 1, 3/2 . . .)-, a basis is given by vectors that, in ket notation, are represented as 
\L,M) -with M = (-L, —L + 1, . . . , L — 1, L)-. One has: 

u € SU{2) R £ b(c n ^J . (2.1) 
The matrix elements of this representation are given by: 

(L, M\RW (u) \L, M') = D L MM , (u) . (2.2) 

These are called Wigner functions [22]. 

The second step is to fix a fiducial state. One can choose the highest weight in the represen- 
tation: \ipo) = \L,L). If the group manifold is parametrised by Euler angles, then u represents 
a point whose "coordinates" range through a £ [0, 4ir) , (3 € [0, n) , 7 £ [0,2-7r) . Fixed the 
fiducial vector, its stability subgroup H^ by the representation is made by elements for 
which (3 = (this condition is seen to be valid independently of N). Two elements u and u' are 
equivalent if u^u' € H^ . It is possible to prove that: 

SU(2)/H i)0 ^S 2 (2.3) 

identifying = (3 and ip = a mod 2ir. Chosen a representative element u in each equivalence 
class of the quotient, the set of coherent states is defined as: 

\9,^,N) = R^(u)\L,L) . (2.4) 

The left hand side ket now explicitly depends on N, the dimension of the space on which the 
representation takes place. Projecting on the basis elements, one has: 

{L,M\0,<p,N) = D L ML {u) , (2.5) 
\0,<P,N) = ( (£+J ff Lm). ) V2 (cos 6/2) L + M (sin 0/2) L ~ M e~* M \L, M) . 

M=-L 

This set of states is nonorthogonal, and overcomplete (dtt = dtp sin d9): 

r "i 2L 

(6',ip',N\6,ip,N) = e ~ iL ^'-^ e*^ - *') cose/2 cos e'/2 + sin 5/2 sin ^/2 

2L + 1 



1 = 

4vr ./ S 2 



[ dn\e,<p,N){0,<p,N\ . (2.6) 
Js 2 



Using this set of vectors it is possible to define a map, from the space of operators on a finite 
dimensional Hilbert space to the space of functions on the sphere S 2 : 



AW G B 



; M N 



#' G T (S 2 ) , 
(9, V ,N\A^\9, V ,N) . 



(2.7) 



This is a way to map every finite rank matrix into a function on a sphere, called the Berezin 
symbol of the given matrix [18]. Among these operators, there are Yj M ^ whose symbols are the 
spherical harmonics, up to order 2L (here J = 0, 1, . . . , 2L and M = — J, . . . , +</): 



{0,<p,N\Y}£\0,<p,N) = Y JM {e,v) , 



(2.8) 



these operators are called fuzzy harmonics. The origin of this name must be traced back to the 
definition, on each finite rank matrix algebra Mjv, of an operator, in terms of the generators 

■Li a • 



^abc 



(2.9) 



representing the Lie algebra of the group SU (2) on the space C 



N 



V 2 A {N) 



M N i — ^ M N , 
Li N \\Li N \AW 



(2.10) 



This operator is called the fuzzy Laplacian. Its spectrum is given by eigenvalues Cj = j (j + 1), 
where j = 0, . . . , 2L, and every eigenvalue has a multiplicity of 2j + 1. The spectrum of the fuzzy 
Laplacian thus coincides, up to order 2L, with the one of its continuum counterpart acting on the 
space of functions on a sphere. The cut-off of this spectrum is of course related to the dimension 
of the rank of the matrix algebra under analysis. The fuzzy harmonics are the eigenstates of 
this operator. 

Fuzzy harmonics are a basis in each space of matrices Mjy. They are orthogonal but non 
normalized and are proportional to the polarization operators Tj^ defined in terms of the 
Clebsh-Gordan coefficients for SU{2) [22]: 
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1 JM 



Sjj'Smm' , 
(-1) M f (iV) 
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the proportionality being: 



y(») 
1 JM 



J! 



(2ijy 



(21 + J + l) 



4vr(2/ - J)! 
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from which: 



Tr 



1 JM I J'M' 



J! 
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47T (21 - J)\ 



Sjj'Smm 1 



Using this basis, an element PW belonging to can be expanded as: 



2L 



j 

E 

J=0M=-J 



jpW y(») 
r JM 1 JM ' 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



with coefficients 



t JM ~ 



Tr 



±l 1 JM 1 JM 



(2.15) 



A Weyl-Wigner map can be denned simply mapping spherical harmonics into fuzzy harmonics: 

(2.16) 



Yjm Yjm(0,v) . 



This map clearly depends on the dimension N of the space on which fuzzy harmonics are realized. 
It can be linearly extended by: 



2L 



f (n) = E E *i N Jt 



JM 



2L +J 



FW(M) = E E **j$YjmM)- ( 2 - 17 ) 



J=0M=-J 



J=0 M=-J 



This is a Weyl-Wigner isomorphism and it can be used to define a fuzzy sphere. Given a function 
on a sphere, if it is square integrable with respect to the standard measure dfi, then it can be 
expanded in the basis of spherical harmonics: 



oo J 



/(0^) = E E fjMY JM (o,<p). 



(2.18) 



J=0M=-J 

Now consider the set of "truncated" functions: 

2L J 



fW(e, ip ) = J2 E fjMY JM {o,<p) 



(2.19) 



J=0 M=-J 



This is a vector space, but it is no more an algebra, with the standard definition of sum and 
pointwise product of two functions, as the product of two spherical harmonics of order say 2L 
has spherical components of order larger than 2L. We have indeed: 

J=\J'+J"\ j 

E I 

J=\J'-J"\ M=-J 



(Yj>m'Yj»m") (0, ip) = E E \l ~ — 4^ (2 J + i) + ~ J "° C J'M' J"M" Yjm {0, if) 



(2.20) 

in terms of Clebsh-Gordan coefficients for SU (2). If these truncated functions are mapped, via 
the Weyl-Wigner procedure, into matrices: 



2L +j 



(2.21) 



J=0 M=-J 



then this set of truncated functions is given the vector space structure of M^r . But M^r is even 
a nonabelian algebra. Invertibility of this association (2.21) enables to define, in the set of the 
truncated functions on the sphere, a non abelian product, isomorphic to that of matrices: 



2L +j 

<7 (JV) )(^) = E E 



J Tr 



J=0 M=-J 



±L 1 JM 1 JM 



(2.22) 
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The Weyl-Wigner map (2.17) has been used to make each set of truncated functions a non 
abelian algebra A^ (S 2 ,*), isomorphic to M^v- The product of two fuzzy harmonics can be 
obtained in term of the 6j-symbols from the product of two polarization operators [7] 



T^jZ" = E (-!) 2L+J V(2 J' + 1) (2 J" + 1) 



J' J" J 
L L L 



n JM rp(N) 

{ ^J'M I J"M" 1 JM ■ 



(2.23) 



These algebras can be seen as formally generated by matrices which are the images of the 
norm 1 vectors in M 3 , i.e. points on a sphere. They are mapped into multiples of the generators 
of the Lie algebra: 



x (l 



■,(N) 



X 



2ie 



VN 2 -1 



(2.24) 



The commutation rules satisfied by generators of the algebras in the sequence A^ (S 2 , *) make 
it intuitively clear that the limit for N — > oo of this sequence is an abelian algebra. This is the 
reason why this sequence is called fuzzy sphere. 



3 The fuzzy disc 

In the previous section, the fuzzy sphere has been introduced via the Weyl-Wigner formalism, 
which makes use of the fact that the 2-dimensional sphere is the coadjoint orbit of the compact 
group SU (2). Compactness of the group brought to a natural definition of a sequence of finite 
dimensional Hilbert spaces, and a natural identification of the dimension of these spaces as a cut- 
off index for a suitable expansion of a generic function on the sphere. Properties of the sphere as 
an orbit of that group made it possible to formalize every point of the sphere as a coherent state, 
defined on each of those finite dimensional spaces carrying the UIRR's. So the map between 
operators (finite rank matrices) and functions on the sphere has been naturally introduced via 
the Berezin procedure, and has been refined stressing the role of the fuzzy Laplacian operator in 
the definition of fuzzy harmonics as a basis in the set Mjv. Moreover, compactness of the group 
played a fundamental role in the rigorous analysis of the convergence performed by Rieffel [6] . 

To approach the problem of defining a fuzzy approximation to the disc, it would be natural 
then to analyse the possibility that the disc were a coadjoint orbit for a Lie group. If this were 
the case, it would be possible to introduce a system of coherent states labelled by its points, 
and some sort of Berezin map [18]. However this approach meets some troubles. The group 
from which it is possible to define a system of coherent states in correspondence with points of 
a disc is SU (1,1), but it is non compact, so its UIRR's are not realized on finite dimensional 
Hilbert spaces. Hence, there is no intrinsic concept of a cut-off index, the would be dimension 
of the fuzzification, in this context. In [14] the problem has been circumvented using coherent 
states for the Heisenberg-Weyl group and the standard isomorphism in terms of Berezin symbols 
between operators and functions on the plane. Then a sequence of projections converging to the 
characteristic function of the disc has been defined. The projections identify a sequence of finite 
rank matrix algebras which can be endowed with an additional structure, a fuzzy Laplacian, 
identifying the underlying geometry. In the commutative limit NO = R 2 this geometry is seen 
to converge to a disc of radius R. 



It is well known that it is possible to define a basis for the space of functions on the disc 
in terms of a suitable set of Bessel functions. Since the Weyl-Wigner map for functions on 
the sphere has been introduced (2.17) mapping spherical harmonics into fuzzy harmonics, we 
pose the question whether or not it is possible to define a system of fuzzy Bessel functions 
in terms of finite rank matrices, and then map again continuum Bessel functions into them. 
Since Bessel functions are defined on the plane, in connection with boundary values problems 
for Laplacian operators, the noncommutative plane which we are going to revue, together with 
fuzzy derivations and a fuzzy Laplacian will turn out to be a natural setting to perform the 
analysis. The final answer to the problem will be given studying the spectral resolution of the 
fuzzy Laplacian. 



3.1 The noncommutative plane as a matrix algebra 

The noncommutative plane can be defined using a Weyl-Wigner map, following again the general 
procedure of Berezin: so the first step will be the definition of a set of generalised coherent states 
for the Heisenberg-Weyl group, since they are labelled by points of a plane. In this space it will 
be assumed a system of coordinates of the kind (x,y) for a point of IR 2 ; 6 will be a parameter 
representing an explicit noncommutativity. 

With this notation, Heisenberg-Weyl group is a manifold IR 3 , whose points are represented 
by a triple (x,y, A), with the composition rule: 

(x, y, A) • (x', y', A') = (x + x', y + y', A + A' - - Q (xy' - x'yfj . (3.25) 

The identity of this group is given by: 

id w = (0,0,0) , (3.26) 

and the inverse of a generic element is: 

A)" 1 = (-x,-i/, -A) . (3.27) 
Considering the complex plane via z = x + iy and ~z = x — iy the group law becomes: 

(z, A) • (z', A') = + z', A + A' + ^ (zz' -zz)^j . (3.28) 

The Hilbert space on which representing this group is again the Fock space T of finite norm 
complex analytical functions in the w variable where the norm is obtained by the scalar product: 

{fl9) = J ^ e_W9 /H5H , (3.29) 
and an orthonormal basis is given by: 

^ {w) = 7W^\- (3 - 30) 

The unitary representation of the Heisenberg-Weyl group is given by operators T (z, A) G 

(f /) H = e ** e -«/2» e W0/ ( w -z) . (3.31) 



Chosen the first basis element ipo (w) as the fiducial state, the procedure already outlined 
gives a coherent state for each point of the complex plane. In the realization of the Fock space 
as a space of complex analytical functions, a coherent state is then: 

\z) -»• i\) z (w) = e -~ zz ' 26 e zw ' e (3.32) 

and, given ip n an element of the basis already considered: 

oo n 

\z) = Y.*-~ ZZ ' 2e -rW^) ■ (3-33) 

n=0 

Such coherent states are non orthogonal, and overcomplete: 

(Z\Z') = e -(\*\ 2 +\zf-2-zz')/0 j 

r <fz 

1 = J^\z)(z\. (3.34) 



On this Hilbert space T, it is possible to introduce a pair of creation-annihilation operators: 

(a/)H = ef- 

aw 

(at/) (w) = wf (w) , (3.35) 

such that: 

a, at] =91, (3.36) 

and: 

a\z) = z\z) , (3.37) 

so that we can recognize the usual coherent states encountered in quantum mechanics textbooks. 
In the chosen orthonormal basis, one has: 



at|y>„) = y/(n + l)6\^ n+1 ) , (3.38) 

and therefore one recognizes the basis as the basis of eigenvectors of the number operator N = 
wa. Note however the unconventional presence 6 in the commutation relation (3.36). 

These relations can be extended. A Berezin symbol can be associated to an operator in the 
Fock space: 

f(z,z) = (z\f\z) . (3.39) 
This can be seen as a Wigner map. It can be inverted: 

f = [ [ flf(z,z)e-fr-i*V 9 ^/ 9 e-& e . (3.40) 
/ nd J ttO 

This quantization map for functions on a plane can be given an interesting form. To start 
with, we can restrict to functions which can be written as Taylor series in z, z: 

oo 

f(z,z)= £ fl7z m z n . (3.41) 

m,n=0 
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An easy calculation shows that this / is the symbol of the operator: 

00 

/= E fln* ]m a n • (3-42) 

m,n=0 

More generally we can consider operators written in a density matrix notation: 

00 

/= E fmn\^m)(^n\ • (3.43) 
m,n=0 

The Berezin symbol of this operator is the function: 

f(z,z) = e-\'V° £ / mB g ^ , (3.44) 
m ,n=0 VmlnW m + n 

where the relation between the Taylor coefficients fmrF and the f mn is 

flk= 2^ fl-qk-q Wq ■ (3-45) 

while the inverse relation is given by: 

min(m,n) p 
fmn V = z2 1 / / x, 7 M/im+n ^ m -P n ~P ' (3.46) 

Equation (3.42) shows that the quantization of a monomial in the variables z, z is an operator 
in a,a^, formally a monomial in these two noncommuting variables, with all terms in a) acting 
at the left side with respect to terms in a. This ordering is related to the presence, in the 
quantization map (3.40), of a specific term that, in the standard terminology used for the Weyl- 
Wigner formalism, is called weight [23]. In fact, restoring real variables (3.40) can be written 
as: 

t f dadb f dxdy , . . j(tx-.») „ a t oa f dadb ~ . ( a2 +>> 2 ) (" &t -" a ) 

/ = / 2^ / ^Jf^y)e-^^^=J ^/(M)e^e^. (3.47) 

with u = (a + z6) /2. The function / is the symplectic Fourier transform of the function /. 
Compared (3.47) with the standard Weyl map [23, 21], the difference is the factor e ( a2+b2 )/ se . 

The invertibility of the Weyl map (on a suitable domain of functions on the plane) enables 
to define a noncommutative product in the space of functions, known as Voros product [24, 25], 
a variant of the more popular Gronewold-Moyal product [26, 27]: 

(f*g)(z,z) = {z\fg\z) . (3.48) 

It is a non local product: 

(/ * g) (z, z) = e~^ e J ti f (z, g (£ z) e^'We*' 9 . (3.49) 
Its asymptotic expansion, on a suitable domain, acquires the form: 

(f*g)(z,z)=fe e ^'g, (3.50) 



and makes it clear that it is a deformation, in 6, of the pointwise commutative product. Since it 
is the translation, in the space of functions, of the product in the space of operators, if symbols 
are expressed in the form (3.44), then the product acquires a matrix form: 

oo 

(/ * 9) m n = Yl fmk9kn ■ (3.51) 
fc=0 

The space of functions on the plane, with the standard definition of sum, and the product 
given by the Voros product (3.48), is a nonabelian algebra, a noncommutative plane. This 
algebra Ag = (^F (M, 2 ^j is isomorphic to an algebra of operators, or, as equation (3.51) 
suggests, to an algebra of infinite dimensional matrices. A more detailed mathematical analysis 
of the Moyal algebra, i.e. of the quantum plane, is given in [28] and references therein. 



3.2 A sequence of non abelian algebras 

A fuzzy space has been presented as a sequence of finite rank matrix algebras converging to an 
algebra of functions. The meaning of this convergence is formalised as compact quantum metric 
spaces. In the case of the fuzzy sphere the rank of the matrices involved is the dimension of 
the Hilbert space on which UIRR's of the group SU (2) are realised. In the approach sketched 
in the last section there is no natural definition of a set of finite dimensional matrix algebras. 
This is a consequence of the fact that the noncommutative plane has been realised via a Berezin 
quantization based on coherent states originated by a group, the Heisenberg-Weyl, which is 
noncompact. 

In this context the strategy to obtain finite dimensional matrix algebras is different. A$ can 
be considered, once a basis in the Hilbert space has been chosen, as a matrix algebra made 
up by formally infinite dimensional matrices. One can define a set of finite dimensional matrix 
algebras simply truncating Ag- The notion of truncation is formalised via the introduction of a 
set of projectors: 

N 

^ (A0 = EKWn| (3.52) 

n=0 

in the space of operators. Their symbols are projectors in the algebra Ag of the noncommutative 
plane, in the sense that they are idempotent functions of order 2 with respect to the Voros 
product (here z = re lip ): 



, N r 2n 
n=0 n=0 



P { e N) {T^) = Y J ^n)^n\z)=e- r2 l 6 Y J 

P e {N) * P W = pf ) . (3.53) 
This finite sum can be performed yielding a rotationally symmetric function: 

in terms of the ratio of an incomplete gamma function by a gamma function [29]. If 9 is kept 
fixed, and nonzero, in the limit for N — > oo the symbol Pg N ^ (r, ip) converges, pointwise, to 



o 



the constant function Pjj N ^ (r,<p) = 1, which can be formally considered as the symbol of the 
identity operator: in this limit one recovers the whole noncommutative plane. 

This situation changes if the limit for N — > oo is performed keeping the product N9 equals 
to a constant, say R 2 . In a pointwise convergence, chosen R 2 = 1: 



1 r < 1 
1/2 r = l 
r > 1 



Jd(r) . (3.55) 



This sequence of projectors converges to a step function in the radial coordinate r, the char- 
acte 
by: 



acteristic function of a disc on the plane. Thus a sequence of subalgebras Ag N ^ can be defined 



4 N) = P r*Ae*Pr ■ (3-56) 

As it has been said, the full algebra Ao is isomorphic to an algebra of operators. What the 
previous relation says is that A^ is isomorphic to Mjv+i, the algebra of (N + 1) rank matrices: 
the important thing is that this isomorphism (this truncation) is obtained via a specific choice of 
a basis in the Fock space T on which coherent states for the Heisenberg-Weyl group are realised. 
The effect of this projection on a generic function is: 



m ^lo VmlnW m + n 

On every subalgebra Ag N ^ , the symbol Pg N ^ (r, tp) is then the identity, because it is the symbol 
of the projector p( N ) = J2n=o IVViX^nl, which is the identity operator in A$ N \ or, equivalently, 
the identity matrix in every Mjv+i- 

Note that the rotation group on the plane, SO (2), acts in a natural way on these subalgebras. 
Its generator is the truncated number operator N( N ) = J2n=o n $IVVi)(^n|- Cutting at a finite 
N the expansion provides an infrared cutoff. This cutoff is "fuzzy" in the sense that functions 
in the subalgebra are still defined outside the disc of radius R, but are exponentially damped. 
In general functions are close to their projected version fi N ^ if they are mostly supported on 
a disc of radius R = VN6, otherwise they are exponentially cut, provided they do not present 
oscillations of too small wavelength (compared to 0). In this case the projected function becomes 
very large on the boundary of the disc. More details and examples are in [14, 15]. 



3.3 Fuzzy derivatives and fuzzy Laplacian 

So far the Weyl-Wigner formalism, and the projection procedure, have provided a way to asso- 
ciate to functions on the plane a sequence of finite dimension {N + 1) x {N + 1) matrices. An 
appropriate choice of 9, the noncommutativity parameter introduced by the quantization map, 
and N, showed that it is possible to obtain a good approximation of a certain class of functions 
supported on a disc. The next step is the analysis of the geometry these algebras can formalize. 
To pursue this task, we need derivations and a Laplacian. The starting point to define the 



i n 



matrix equivalent of the derivations is: 

dzf 

dzf 



/,« f \z) , 



a,f \z) . 



(3.58) 



This relation is exact in the full algebra Ag. Given an operator /, the derivatives of the symbol 
/ (z, z) are related to the symbol of the commutator of / with the creation and annihilation 
operators. A fuzzified version, namely a truncated version, of these operations, can be defined 



as 1 : 



i 



9 >' 



dzf. 



(JV) 



- - l MH N) \H N) !H N \a 



H N) \ Z ) 
H N) \z) • 



It is important to note that this is really a derivation on each A g N \ that is a linear operation 
from A g N ^ to itself, satisfying the Leibnitz rule. The way this is achieved in (3.59) goes through 
the following steps. We first associate to fi N ^ a finite rank matrix fjj N ^ which is indeed regarded 
as embedded in an infinite dimensional one, with all zeroes but (N + 1) x (N + 1) elements in 
the upper left corner: f@ N ^ = Pg N ^ fPg N \ We then perform the commutator with fi, fit realised 
as infinite dimensional matrices. This operation yields an infinite dimensional matrix with all 
zeroes but (N + 2) x (N + 2) elements in the upper left corner, so that we have to project back 
to Mat + i to obtain a properly defined derivation. Finally, we consider the symbol. 

As an example let us calculate the derivatives of the fuzzified coordinate functions z and z. 
From the definition, z, z are respectively the symbol of the annihilation and creation operators, 
which can be written as: 



(3.59) 



s=0 

oo 

= V(k + i)e\4> k+1 )(i> k \ 



(3.60) 



Projection into Mat + i gives: 



fc=0 



N-l 



(3.61) 



s=0 



According to (3.59), to perform the derivative with respect to z, we have to consider: 

N-l 



.(AT) .+ 



\ipo){M + Yl " N\i(> N ){il> N \ 

s=l 

N-l 

Y\4>s)(4>s\-NW N )(^ 



s=0 
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s=0 



(3.62) 



1 This a slightly modified version of the derivations in [14, 15] 
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The first term of the sum is the projector onto the first N + 1 basis elements, whose symbol is 
the identity in A$ N \ It is worth noting that this commutator has no terms 'outside' the space 
we are considering, namely there are no components on density matrices of order greater than 
N + 1: this means that in this case there is no need to project it on Mjv+i- The symbol yields 
then: 



d z (z. 
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^ —ttt— r 2N N N 

ls=0 Sl 



(3.63) 



In the limit of N — > oo the first term is the characteristic function for the disc, while the second 
converges to a factor tt5 (r — 1). This factor is a radial 5 selecting the value r = 1 with respect 
to the Lebesgue measure on the plane: 



lim d z (z { N) ) = Id (r) - ttS (r - 1) . 



N 

To calculate the derivative of z with respect to z one needs to consider: 

JV-1 



(3.64) 



4 ( ) = E VW+voWk+i)^ 
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and the commutator: 



-0y/N(N+l)\1> N+1 ){lJ> N \ . 



(3.65) 



Unlike the previous example this operator must be projected back to the algebra Mjv+i, and 
upon considering the symbol we finally have: 



)-»■ 



(3.66) 



Let us come to the definition of the Laplacian operator. In the spirit of noncommutative 
geometry it is this additional structure which carries the information about the geometry of 
the space underlying Ae- In facts, we can say we have succeeded in fuzzifying the disc only if 
we have been able to define a fuzzy Laplacian whose spectrum approaches that of the ordinary 
Laplacian on the disc when N — ► oo. 

Starting from the exact expressions: 



V 2 f(z,z)=4d- z d z f = -,(z\ a, /,ot \z) 



it is possible to define, in each Ag N ^: 
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(3.68) 



The image of the element of the truncated algebra: 

N 

f(N) 



a,b=0 



1 



is then: 
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-EE^( s + x ) i^-x^i - E /o.«+i ( s + x ) i^oX^+ii + 

s=0 6=0 s=0 
N-l N-l 

+ E E + + 

s=0 6=0 
JV-1 AT-1 
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(3.69) 



The eigenvalues of this Laplacian can be numerically calculated. They are seen to converge to 
the spectrum of the continuum one for functions on a disc, with boundary conditions on the edge 
of the disc of Dirichlet homogeneous kind [30]. All eigenvalues are negative, and their modules 
A solve the implicit equation: 

Jn = , (3.70) 

where n is the order of the Bessel functions. In particular, those related to Jo are simply 
degenerate, the others are doubly degenerate: so there is a sequence of eigenvalues labelled by 
\ n ^k where n is the order of the Bessel function and k indicates that it is the k th zero of the 
function. The eigenfunctions of this continuum operator are: 



with n integer number and |n| its absolute value. This is a way to write the eigenfunctions in a 
compact form, taking into account the degeneracy of eigenvalues for \n\ > 1. 

The spectrum of the fuzzy Laplacian is in good agreement with the spectrum of the con- 
tinuum case, even for low values N of the dimension of truncation, as can be seen in figure 1. 
It is interesting to note that with this definition of Laplacian (3.68) we correctly reproduce the 
pattern of non degenerate and double degenerate eigenvalues. The difference with the case of the 
spectrum of the fuzzy Laplacian for the fuzzy sphere is that now the "fuzzy spectrum" is both 
a cut-off and an approximation of the continuum spectrum. It is a cut-off because, of course, it 
is a finite rank operator. It is an approximation because it has been defined using a formalism 
whose building blocks are related to a noncompact group, namely the Heisenberg-Weyl, whose 
generators have no finite dimensional realization. 



4 Fuzzy Bessel functions 

In this section we introduce a basis on the fuzzy disc. As we have already noticed in the 
introduction, we cannot use representation theory, which in the case of the fuzzy sphere yields 
the fuzzy harmonics, because there is no compact group whose homogeneous space is the disc. 
But we have a well defined Laplacian for each Ag N \ Therefore we may answer the problem 
looking at the eigenfunctions of this fuzzy Laplacian. 
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Figure 1: Comparison of the first eigenvalues of the fuzzy Laplacian (circles) with those of the 
continuum Laplacian (crosses) on the domain of functions with Dirichlet homogeneous boundary 
conditions. The orders of truncation are N = 10, 20, 30. 



4.1 Spectral properties of the Laplacian 

The fuzzy Laplacians (3.68) are a family of operators mapping each algebra into itself with 
the explicit action (3.69). Considering Ag N ^ as a (N + l) 2 dimensional vector space, a density 
matrix |V'a)(V ; fe| is a basis element for it, with a, b G (0, . . . , TV). Equation (3.69) can be specified: 

V^dVaX^I) = 4ArP e {A0 (v^a-iX^-il + (4.1) 
-(a + 6 + 1) m(ipb\ + V(a + l)(b+l)\4, a+1 )(4; b+1 \)Pi N) . 

The action of the projector operator is to cut away some of the components of the right side 
image for a = 0,N and for b = 0, N. It is important to note that, chosen an integer n = a — b, 
the fuzzy Laplacian maps basis elements with the same n onto themselves, as n ranges through 
(— N, . . . , +N). In the matrix representation of A$ N \ the number n fixes the diagonal to which 
the element \ip a }(^b\ belongs. The equation (4.2) then shows that the fuzzy Laplacian maps 
each diagonal subspace of M^+i into itself. This property simplifies the study of the spectral 
properties of the fuzzy Laplacian. The analysis of these spectral properties will be performed 
restricting the operator to each of its stability subspaces. 

To understand the meaning of this integer n, go back to the eigenfunctions of the continuum 
Laplacian on a disc with Dirichlet homogeneous boundary conditions (3.71). Eigenfunctions 
&n t k are defined on the real plane, so they can be mapped into operators on the whole Hilbert 
space via the Weyl map introduced before (3.40). The continuum eigenfunctions (3.71) with a 
fixed nonnegative n, such that < n < N can be mapped into an operator: 

(¥)" /2 §(-¥)* jttW 6, ' s " + "- <"> 

In the density matrix notation, it acquires the form: 

s\ (s + n)\ (j — s — n)\ 
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This element can be fuzzified, so that now the parameter 6 is constrained as ON = 1: 
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n.fc 
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i/2 /V- 



y/a\{a + n)\ 



a=0 
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,s=0 
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n,k 



AN J s\ (s + n)\ (a - s)\ 
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(4.4) 



Eigenfunctions with negative n are just the complex conjugate of those with positive n: ($_ ni fe)* = 
$ n> fe. The fuzzification procedure gives a matrix which can be obtained by by an hermitian 

conjugation: = ^j^- This analysis clarifies that, in the fuzzy approximation, continuum 

eigenfunctions are represented by matrices belonging to the n th diagonal subspace of M/y+i- 
Moreover, in this approach, admissible n are constrained by the dimension of the fuzzification: 
n < N. In all these fuzzy elements, A„£ (as well as Aofc) can be seen no longer as eigenvalues 
of the Laplacian operator on the disc with Dirichlet homogeneous boundary conditions, but, at 
this stage, simply as parameters. 

We can now perform a complete spectral analysis of the fuzzy Laplacian. This will also show 
the exact meaning of the X n> k parameters. Consider first the n = subspace, whose dimension 
is N + 1, and consider the elements with k = 1, . . . N + 1. Act on these elements with the 



fuzzy Laplacian. From equation (4.2) we obtain: 
V^WaWal) = 4iV (a|^o-i>W«-i| - (2a + 1) |^)(Va| + (a + 1) IVo+lWa+ll) 
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(/V) 



(4.5) 

^rffc ) on tne basis density matrices for the specific subspace. 



one can reproject the image V?^ 
On the first basis element, equation(4.5) becomes: 

<Y>o| V 2 {N) (m&a\) I'M = AN(ad , a -i - (2a + 1) S 0>a ) 
Considering the explicit form of (4.4) we obtain, for n = 0: 



(4.6) 
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where: 



\k 



and upon applying the Laplacian: 
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= -A ,fc = -A ,fe (^o|*o,fc IV'o) 

(TV) 



4iV - 
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(4.9) 



Hence the fuzzy Laplacian maps the element <3?q k into another diagonal element whose coef- 
ficient, respect to \'4>o)(4>o\, is just a multiplication by — Ao,fc of the coefficient of respect 
to the same basis element. Proceeding along this way we project the element (&l^ k ^ 



on 



1, 



,N 



\fpb){fpb\, with b 



h (b,N) 
°0,k 



) + 6(<-^'-<f). (4.10) 
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To study the r.h.s. of this relation, we need to calculate the difference: 
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_ a (6,JV) 
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hence the equation (4.10) becomes: 



(H v (;V) v ■ 



(TV) 



Mb) 




(b-l,N) 
0,k 
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(4.12) 



where to reach the expression on the last line we have recollected equal powers of the A and made 
some simplifications. This relation is similar to (4.9). Whatever Ao,fc, the fuzzy Laplacian acts 
on the first N components of the fuzzy elements iff simply as a multiplication by — Ao,fc- To 
conclude that the elements iff are eigenstates of the fuzzy Laplacian we need to study the last 

component projection of (iff), namer y that, on the basis element \iPn)(iPn\- Moreover, 

we have not yet obtained any conditions on the possible eigenvalues of the fuzzy Laplacian. 

Let us project: 



<^v| Vfo (iff) |^> = AN (N*$-W - (2N + 1) 



and impose: 

<lM Vfo (iff) \^n) 
This can be cast in the form: 
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Recollecting terms with equal powers of A, the equation can be written as: 
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(4.15) 



(4.16) 



The meaning of this calculation is clear: relation (4.14) is not identically satisfied by all value of 
Ao,fc. So it becomes an equation, giving the admissible eigenvalues Aff of the fuzzy Laplacian 

on the subspace of diagonal fuzzy elements in Ag 
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if and only if: 
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(4.18) 



We now show that the solutions of this last equation are actually the eigenvalues of the fuzzy 
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subspace, for $(^0 = J2a=o 0^ IV'aKV'al we have: 



Laplacian (in the n = case). From the action of V^n on a generic element of this diagonal 
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the eigenvalue equation (4.17) can be written in components as: 



A, 



(JV) 
0,fc 



4iV 



(a + 1) 



(TV) 
a+1 



+ 



^-(2^ + 1)^ 



(TV) 
■Q,fc 

4iV~ 

(JV)' 

ITT 



A, 



A 



5 JV 



if a / 0,iV 



(4.20) 



This set of equations can be given a matrix form. If we set xjv = Aq^/4AT, then the eigenvalues 
are related to the roots of the characteristic polynomial Pn (xjv), given by the determinant of 
the matrix 



P/v = det 



/ x N - 1 1 
1 xn — 3 
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N x N -(2N + l) ) 



The eigenvalue equation for the fuzzy Laplacian in the diagonal matrices subspace is: 
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(4.22) 



Matrices of the kind of (4.21) ara called Jacobi matrix [31], and their characteristic polynomial 
is given inductively: 



P/V (x N ) = (X N ~ (2N + 1)) Pv-l (X N ) ~ N 2 Pjv-2 (x N ) ■ 



(4.23) 



Equation (4.22) is equivalent to (4.18). This result, proven in the appendix, closes the analysis 
on the spectral properties for the n = case: eigenstates are given by $qj?' wm l e eigenvalues 
are given by the roots of (4.22). 
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We approach the problem for the action of the fuzzy Laplacian in the nondiagonal subspaces 
of Ag N \ defined by the integer n, following a similar strategy. The dimension of a generic 
subspace is N + 1 — n, the basis elements are \ip a ) {4> a +n\ with a = . . . , N — n so the fuzzy 



Laplacian can be seen as an operator acting on C 



N+l-n 



Going through the path described for 



the n = case, we prove that the eigenvalue problem in this subspace is solved by the fuzzy 



elements (4.4) with k = 1, . . . TV + 1 
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(4.25) 



This equation generalises the equation (4.16), reducing to that for n = 0. 

The proof will closely follow the path outlined in the diagonal case. First project the image 
element V 2 ^ (&^^) on the first basis element |V'o)(V'n|- Following the general case (4.2), one 
can see that: 



V*N) (|^o)<^a+n|) l^n) = 4A [S^yfr+T - {n + 1) 6 0>a ] 



After some algebra using the definition (4.4): 
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where the value \ n ^k is again just a parameter. 
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0,N — n. Again, after some straightforward algebra: 
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This r.h.s. defined to be: 
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This is an equation for X Ht k, that can be cast exactly in the form of (4.25): 
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So the eigenvalues of the fuzzy Laplacian in this stability subspace are given by the solution 
of: 
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(4.32) 



It is again possible to prove that this equation is exactly equivalent to the secular equation 
coming from the matrix representation of the action of the fuzzy Laplacian on each of this 
N — n + 1 dimensional subspaces. The analogue of the matrix (4.21) is, in this case: 
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N — 



( x N -Jn+l ) y / (n + 1) 
y/(n + 1) x N -(n + 3) 
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y/N (N — n) x N - {2N-n+l) J 



(4.33) 

The characteristic polynomial is of degree is N + 1 — n and a recursive relation holds, analogous 
to (4.23): 



W ~ (x N - (2N - n + 1)) P£ (x N ) -N(N-n) p£ ( Sjv ) 



In the appendix there is the proof of the equivalence of the relation 







(4.34) 



(4.35) 



with (4.32). Setting n = the diagonal case is recovered. 

The analysis for the case n < is straightforward. The eigenvalue equation is the same, and 
this shows why there is a set of doubly degenerate eigenvalues. The eigenstates are obtained by 
those for positive n by complex conjugation. This completes the spectral analysis of the fuzzy 
laplacian. Since we have proved that the eigenstates of this fuzzy Laplacian are the fuzzifled 
version of the continuum eigenfunctions, we call fuzzy Bessels. 



4.2 Comparison of Bessel Functions 

We now compare the behaviour of the symbols of the fuzzy Bessel with their ordinary counter- 
parts. From (4.4) the symbol of a fuzzy Bessel is: 
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(4.36) 



The integer n appears as a phase modulating factor for the variable <p. This would be the 
expansion of the corresponding Bessel function, where it not for the truncation in the sum, and 
the fact that the parameter Ao,fc has become the eigenvalue of the fuzzy Laplacian, i.e. a solution 
of (4.22). For the n = the expression can be simplified: 




(4.37) 
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Where C 



(a) 



(N) 



(JV) 



-jj^- j is the a Laguerre polynomial in the variable J ■ We can plot the 

diagonal fuzzy elements. Fig. 2 shows that the zero order fuzzy Bessel state converges to the 
continuum eigenfunctions $0,1 ( r > f° r values of r inside the disc of radius 1, while it converges 
to zero outside the disc. This behaviour is seen to be valid also for eigenstates of different 






Figure 2: Comparison of the radial shape for the symbol <3?q 1 (r, ip) (continuum line), the symbol 
of the eigenmatrix of the fuzzy Laplacian for N = 10, 20, 30, with $0,1 ( r i <p)- 

eigenvalues. The plots for the symbols are ™ fig ure 3, those for ^q^? in figure 4. It 

n=0,k-2,N=10 n=0,k=2,N=20 n=0,k-2,N=30 
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Figure 3: Comparison of the radial shape for the symbol & 2 (r) (continuum line) with <3?o,2 (r) 
for N = 10, 20, 30. 

is interesting to analyse the fuzzification of $0,10 ( r )- The fuzzy symbol is ^oxq- ^ or ^ = 
10,20,25 it is plotted in figure (5). It is evident that the fuzzy eigenfunction reproduces the 
continuum eigenfunction for values of r close to the centre of the disc, but not on the edge, 
where a huge bump appears. In [14] it has been explained that the presence of the bump 
on the edge of the disc, in the fuzzification of a function defined on the plane, is related to 
the fact that this function has oscillations of too small wavelenght compared to 9. This is 
manifestation of the infrared-ultraviolet mixing characteristic of noncommutative theories. In 
the case of <3?o,io (r) = Jo (y/^o,iof) , one can immediately see that the oscillation wavelenght 
of the continuum eigenfunction is given by p\ ~ 1/k ~ 1/10. In these plots, it is assumed 
6 = 1/N (the fuzzy disc truncation), so 9 and p\ are of compatible magnitude. In the fuzzy disc 
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Figure 4: Comparison of the radial shape for the symbol <&03 ( r ) (continuum line) with <3?o,3 (r) 
for N = 10, 20, 30. 
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Figure 5: Comparison of the radial shape for the ^J/io symbol (continuum line)for N = 10, 20, 25. 
The bump goes out of the picture. As N increases, it becomes narrower and narrower. 
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limit, N — ► oo, so 6 is infinitesimal. The bump disappears, as it is shown in the other plots 
(figure 6). 
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Figure 6: Comparison of the radial shape for the $9^10 s D m bol for N = 30, 35, 40. 



The non radial functions follow a similar pattern. Their phases are exactly as the ones of 
their continuum counterparts, while the radial parts are similar. A first plot is in Fig. 7, a 
second one in Fig. 8, where the fuzzification procedure gives again a bump, for small values of 
N. This bump is seen to disappear in the fuzzy disc limit. 






Figure 7: Comparison of the radial shape of the ^ (r, ip) symbol (continuum line) for N 
10,20,30. 



The fuzzy Bessels, being eigenfunctions, provide an efficient way to calculate a fuzzy Green 
function, which we calculated numerically in [14]. Since the fuzzy Laplacian is represented 
by an hermitian finite dimensional matrix, its inverse can be written in terms of a spectral 
decomposition. If a matrix, say M, is hermitian, then its components satisfy the condition 
M a i, = M£ a in terms of transposition and complex conjugation. The eigenvalue problem gives a 
number of real eigenvalues equal to the dimension of the space on which the matrix acts: 



E^f = A (fc M fe) 



(4.38) 



Figure 8: Comparison of the radial shape of the 3>g5 (r, ip) symbol (continuum line) for N 
10,20,30. 



here indicates the a th components of the eigenvector relative to the eigenvalue A^. The 
inverse, if it exists, of the matrix M is a matrix G whose components can be written as: 

G sq = 4 k) v { q k) /\ {k) ■ (4.39) 
k 

If we specialise to the problem under analysis, the notion of eigenvector of components vi N ^ 
with eigenvalue A^ goes into that of fuzzy Bessel matrix, from which it is immediate to 
obtain the symbols. The fuzzy Green function becomes: 



+N N+1 -\ n \ (N) , ,y $ (JV) ( s 
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n=-N k=l 



Since the fuzzy Bessels play a role similar to fuzzy harmonics for the fuzzy sphere algebra, 
we can now make describe the process of approximating the algebra of functions on a disc with 
matrices more precise. In complete analogy with (2.18) and (2.19), if / is square integrable with 
respect to the standard measure on the disc dQ = rdrdcj), it can be expanded in terms of Bessel 
functions: 

+oo oo 

f(r,<p)= E E f^ eimP J H (VW") (441) 

Tl— OO fc = l 

and it is possible to truncate: 

+ N N+l-\n\ +JV N+l-\n\ 

f N Hr, ( p)= E fnk^J\n\ (yfiw*r) = E E /n* *n,fc (r, *>) . (4.42) 

n=-N k=l n=-N k=l 

This set of functions is a vector space, but it is no more an algebra, with the standard definition 
of sum and pointwise product, as the product of two truncated functions will get out of the 
algebra. The mapping from truncated functions into finite rank matrices: 

+N AT+l-H 

/W ^ f y) = ^ ^ fnk$ W (4.43) 

n=—N fc=l 



endows the set of functions with a noncommutative product, inherited from the matrix product, 
which makes it into a non abelian algebra. The formal limit TV — > oo with the constraint NO = 1 
is the abelian algebra of functions on the disc. The sequence of nonabelian algebras is 
what we call the fuzzy disc. 



A Weyl-Wigner Formalism 

In this appendix we give a very short introduction to the Weyl-Wigner formalism used to define 
maps between functions on a "classical" plane and a set of operators on a Hilbert space [23]. 

Given a symplectic real, finite dimensional, vector space (L,uj), where the symplectic form 
to is translationally invariant, i.e. with constant coefficients, the set of Displacement operators 
define a Weyl system, a unitary representation of canonical commutation relations in the expo- 
nentiated unitary version [32] : 

D : L ^U{H) 
D{z + u) = e iuj{z ' u)/2e D(z)D(u) . (A.l) 

Here is a parameter: the Displacement operators can be seen to set a unitary projective repre- 
sentation of the translation group, where the phase factor is related to the sympletic structure 
defined on the vector space L. These operators also define a representation of the Heisenberg- 
Weyl group [33]. On every one dimensional subspace of L (a is a real scalar), via the Stone- von 
Neumann theorem: 

D {az) = e ia6(z)e (A.2) 

and these generators satisfy: 

G (z) , G («)] = i0u (z, u) 1 . (A.3) 

This shows that a Weyl system is a way to formalise a set of canonical commutation relations for 
quantum observables, recovered as hermitian generators of a unitary ray representation. In this 
perspective acquires the role of a noncommutativity parameter. In the standard approach, 
the symplectic vector space is seen as a phase space for the classical dynamics of point particle: 
generators of a Weyl system represent the position and momentum observable for a quantum 
dynamics of point particle. 

Using the definition of a Weyl system, it is possible to define a Weyl map, that is a map 
from functions on the vector space (Llo) (dim L = 2n) to Op(Tt), the set of operators on the 
Hilbert space on which Displacements are realised: 

/ = / J^ff^)D{w) , (A.4) 

where: 

/H = / w&r f {z) e ~ Mz,w)/e (A- 5 ) 

is the symplectic Fourier transform [33] of the function /. The Weyl map can be seen as a 
quantization map for classical observables defined on a phase space: it can be inverted by the 
Wigner map: 

f(w) = Tr \fD^(w)] . (A.6) 
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So the Weyl-Wigner maps define a bijection between a set of classical observables and a set of 
quantum observables. This is the building block of a formalism suited to study both the problem 
of quantizing a classical system and performing a classical limit of a quantum system. 

The standard Weyl-Wigner maps can be modified by the introduction of a term, called 
weight: 



fx 



I 



dw 



(2tt0)' 



f(w)X(w)D(w) 



f(w) = \(w)- l Tr f x DHw) 



(A.7) 



This weight can be proved to be related to a more general class of ray representation of the 
translation group, and is seen to be responsible both for a specific "ordering" in the quantization 
procedure and a definition of a specific domain of mathematical applicability of the bijection, 
as stressed in the main text. 



B Generalised coherent states 

In this paper the concept of generalised coherent states has been extensively used. The aim of 
this appendix is to briefly recollect the main definitions and results, to fix notations and facilitate 
the reading. The main reference for coherent states is [34]. 

Consider a Lie group G, and U (g) a unitary irreducible representation of this group on a 
Hilbert space TL. Chosen a fiducial vector |V>o) in Ti, one obtains a set of vectors for each element 
of the group, acting on it with U (g): 

ty g ) = U{g)\iM . (B.l) 

Two such vectors are considered equivalent if they correspond, quantum-mechanically, to the 
same state, i.e. if they differ by a phase. So \ip g ) ~ if \ip g ) = e 1 ^ 9 ' 9 '^ \ip g '). This condition 
is equivalent to U (V'o) = e 1 ^ 9,9 '^ \ipo)- If H is the subgroup of G whose elements are 

represented, by U, as operators whose action on the fiducial vector is just a multiplication by a 
phase, then the equivalence relation is among points of G, and the quotient is the space G/H. If 
H is maximal, then it is called isotropy subgroup for the state | V'o) • Choosing a representative 
g(x) in each equivalence class x € X = G/H (which is a cross section of the fiber bundle 
G with base X) defines a set of vectors on H, depending, clearly, on G and \ipo). This set of 
states is called a system of coherent states for G. The state corresponding to the vector \x) 
may be considered as the range of a rank one projector in Ti. Thus, the system of generalised 
coherent states determines a set of one dimensional subspaces in Ti., parametrised by points of 
the homogeneous space X = G/H. An evolution of this analysis drives naturally to the issue of 
overcompleteness for the system of coherent states, mentioned in (3.34), (2.6). 



C Explicit Calculations for the Eigenvalues 

The aim of this appendix is to show an explicit proof of the results claimed in the main text 
about the eigenvalue problem for the fuzzy Laplacian. In particular, we want to prove that the 
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equation (4.32) is equivalent to the equation (4.35). From this equivalence, it will be proved 
also the equivalence between equations (4.18) and equation (4.22), as they are a special case of 
the first two for n = 0. 

The polynomial E$ is given, identifying X { ^/4N = x N , by (4.25): 

^ (^)- 2. fc!(n + fc)!(JV-n-A: + l)! ^ "°- (C - 2) 

The characteristic polynomial is given by the determinant of the matrix (4.33) given re- 
cursively by (4.34). To prove the equality of the two polynomials, the first step is to prove that 
E^ satisfies the same recursive relation that pffl does. 



Considering the quantity: 



(x N — (2N — n + 1)) (x N ) - N(N-n) E^_ 2 (x N ) = 

E'~ n ( u k+N-n N\(N-n)\ fc+1 

( q 1 ' k\{n + k)\{N -n-k)\ { N) + 

-g(-i)^^-» +1 ) M(B ^-:!_ t) , (^ + 

Where the coefficient of the zeroth order term is: 

( _ 1} iV-n + l m {m _ n + 1} _ { _ 1)N -n + l M {N _ \ = ( _ 1)JV -n+l (^+1)! (c _ 4) 

(n) 

It coincides with the coefficient of the zeroth order term of the polynomial E K N ' . 

Let us consider, in the expression (C.3), the coefficient of the q order term, with 1 < q < 
N — n — 1. It is equal to: 

(-l)^-i AT! (N - n)\ ( (a _ 1)!(ra+a _ 1 1 )!(jV _„_ a+1)! + (C5) 

2N-n+l 1 \ („\q _ ,_ U N-n+a+l (N+l)\(N-n+l)\ , sn 

a\(n+a)\(N-n-a)\ a\(n+a)\(N-n-a-l)\ J [ - XN ' ~ v L > a\{a+n)\(N-n-a+l)\ ^ N > 

This coefficient is equal to the coefficient of the q th order term in E^\ The equality of the 
coefficients, of the remaining order terms, of the polynomial (C.3) with those of the polynomial 

N 



E*f* ( x n) can be easily checked. 



At this stage, we have proved that the polynomial E^ (xn) satisfies the same recursive 
relation that the polynomial pffl (ajjv) does satisfy. The recursive relation they satisfy is such 
that the N th term of the sequence E^\ for each fixed n, depends on both the (N — l) th and 
(N — 2) terms of the sequence. To prove the complete equivalence of these two polynomials, 
we thus need to prove that they explicitly coincide at a pair of consecutive steps. 

In the main text it has been stressed that there is a constraint for allowed positive n, namely 
an upper bound: n < N + 1. This means that, for a fixed n, that is for a fixed subspace in 



the fuzzy algebra, N > n+1. For N = n — 1 both polynomials E^ l } 1 (xjv) and (xjv) are 
trivial, so we study at first the case of N = n. We have: 



*•> - E (-^' s^^gfa - * m - <" + 1 > « c - 6 > 

while the matrix representing the action of the fuzzy Laplacian is one-dimensional, so that: 

p( N ) = XN - (2n + 1) . (C.7) 
The second case is for N = n + 1. We have: 

' k\(n + k)\(2-k)\ { N) 
= (x N ) 2 -2{n + 2)x N + (n+l)(n+2) . (C.8) 



The matrix representing the action of the fuzzy Laplacian is now two-dimensional, so that the 
characteristic polynomial is: 

pM = detf ^ , ) = 

= (x N f -2(n + 2)x N + (n + l)(n + 2) . (C.9) 
This proves the equivalence of the two polynomials. 
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